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' ABSTRACT 
ON ' 

0^ ■ As was shown by Leutwyler and Smilga, the fact that chiral symmetry is broken 

and the existence of a effective finite volume partition function leads to an infinite 
number of sum rules for the eigenvalues of the Dirac operator in QCD. In this paper 
we argue these constraints, together with universality arguments from quantum 
chaos and universal conductance fluctuations, completely determine its spectrum 
near zero virtuality. As in the classical random matrix ensembles, we find three 
universality classes, depedending on whether the color representation of the gauge 
group is pseudo-real, complex or real. They correspond to SU (2) with fundamental 
J> i fermions, SU(N C ), N c > 3, with fundamental fermions, and SU(N C ), N c > 3, 

' with adjoint fermions, respectively. 

O ' 
O 

■ 1. Introduction 

o 

^ ■ It is widely believed that QCD is the theory of strong interactions which describes 

the hadronic world as we know it. On the other hand, at low energies, the same physics 
can be described much more economically in terms of a partition function involving 
effective degrees of freedom. This raises the question to which extent the low energy 
effective theory put constraints on the underlying microscopic theory. 

This question can be answered if we restrict ourselves to a space time volume 
V4 = L 4 with L <C m' 1 . In that case only the static modes of the Goldstone fields 
contribute to the partition function. If, at the same time, L 3> Aq^ d , the low-energy 
effective paritition function coincides with the QCD partition functional. Its mass 
dependence is given byi 

Z eff (m,6)= [ DU exp(\/ 4 SReTr me i6/N ^U) (1) 

JUdG/H 

where the integral is over the coset G/H. For the standard scheme of chiral symme- 
try breaking we have G/H = SU(Nf) x SU(Nf)/ SU(Nf). In agreementB with the 
Vafa-Witten i theorem, this corresponds to a condensate £ that is flavor symmetric 
(see also refsfli). For simplicity we take the mass matrix m diagonal in this paper. 
As dictated by QCD Ward identities, the masses always occur in the combination 
me %6 ' N f , where 9 is the vacuum angle. 

The mass dependence of the QCD partition function is given by 

2 QCO Kfl) = E^(Il m / II ^l + m}))u, (2) 

v f A„>0 



where the average (■ • •)„ is over all gauge configurations with topological charge v 
weighted by the gluonic action. 

As was argued correctly by Leutwyler and SmilgaEJ, the the finite volume partition 
function, imposes constraints on the eigenvalues of the Dirac operator. More precisely, 
sum rules are obtained by equating the coefficients of the expansion in powers of the 
quark masses of both sides of the equation 

Z$ CD (m) Zf(m 



efit 

(3) 




Z, QCD (0) Zf (0) ' 

where Z u (m) is the z/'th Fourier coefficient of Z(m, 9). Since this leads to an infinite 
set of sum rules, the following question should be raised: to what extent is the 
spectrum of the Dirac operator determined by the finite volume partition function? 
Because we study a region where only the static modes contribute to the finite volume 
partition function, we expect that it only determines the part of the spectrum near 
zero virtuality . To make this statement more precise let us consider the simplest 
sum rule obtained by equating the 0{m 2 ) terms in @ 

Y 2 

= (4) 

„,„ m+»y 

where the sum is over the positive nonzero eigenvalues of the Dirac operator. If the 
spectral density is defined by 

p(X) = (£'8(\ - A n )> (5) 
we can rewrite the l.h.s of (|J) as 

1 ,v^/ 1 . f , 1 . z s 1 

where we have introduced a new integration variable by 

z = XV 4 . (7) 
This leads to the definition of the microscopic limit of the spectral density! 

It is clear that information of the static finite volume partition function is contained 
in this limiting function of the spectral density. Since Z cS is based on symmetries 
and chiral symmetry breaking only we arrive at the conjecture that the microscopic 
spectral density is completely determined by the symmetries of the QCD partition 
function. The same applies to microscopic correlation functions of the spectral density 
which appear in higher order sum rules. 
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It is our program to contruct a spectral density that contains nothing else but 
the symmetries of the QCD partition function as input. Such question has been 
phrased in general terms by BalianB, who introduced the maximal entropy approach 
to obtain the probability distribution of the ensemble that decribes a system that 
obeys a set of well defined constraints. This approachhas been applied successfully 
: statistical S"-matrix fluctuations™^" and universal conductance 
For example, in the case of the classical random matrix ensembles, 
a Gaussian probabilty distribution of the matrix elements is obtained in this way. 

Along these lines we can construct a random matrix model with the symmetries of 
the QCD partition function as input. As in the classical random matrix ensembleslij, 
we find three different universality classes^!: 1. The chiral orthogonal ensemble (ch- 
GOE) for 577(2) with fundamental fermions when the Dirac operator is real. The 
chiral unitary ensemble (chGUE) for SU(N C ), N c > 3, with fundamental fermions 
(in this case the Dirac operator is complex). 3. The chiral symplectic ensemble 
(chGSE) for SU(N C ), N c > 2, with adjoint fermions (in this case the Dirac operator 
can be regrouped into real quaternions). In all three models the scale is set by the 
chiral condensate, which, via the Banks-Casher formula@, is related to the average 
level density. Further they depend on the number of flavors and the total toplogical 
charge. Precursers of the chGUE discussed in this paper can be found in refs.0@. 

Of course, it makes only sense to use random matrix theory for observables that 
are not sensitive to the detailed dynamics of the system. From studies of rmdearlevel 
Ji§cj. universal conductance fluctuationscIrB3. Quantum chaos EilESE 



resonancesaH, universal conductance fluctuations^^, quantum chaos Ea'Ea'EltBEZl and 
general argumentsHa we know that correlations between eigenvalues on the scale of a 
finite number level spacings are such observables. One of the best known examples is 
the nearest neighbor spacing distribution. Because the micrsocopic spectral density 
(|S|) describes correlations on the scale of a finite number of levels away from zero, 
we expect that this quantity is universal. It is our conjecture that the exact QCD 
microscopic spectral density belongs to this universality class. 

What are the arguments in favor of universality? In the ideal case one would 
like to analyze the spectrum of the QCD-Dirac operator as obtained in numerical 
lattice calculations. However, technically it is not possible to work close enough to 
the chiral limit, and the best one can do in this direction is to work with cooled 
configurations which are represented by a partition function of a liquid of instantons 
and anti-instantonsll (see however^). In this paper we summarize the main numerical 
results (see section 6). Important support in favor universality comes from a different 
branch of physics, namely from the theory of universal conductance fluctuations. In 
that context, the microscopic spectral density of the eigenvalues of the transmission 
matrix, was calculated for the Hofstadter0 model, and, to a high degree of accuracy, 
it agrees with the random matrix prediction^. The Hofstadter model, which has the 
symmetries of the chGUE for zero flavors, is therefore in the same universality class 
as the corresponding random matrix model. Because, this model differs in all other 
respects from the random matrix model, this result illustrates the size of the basin of 
attraction of the random matrix model and greatly increases the hope that the exact 
QCD microscopic spectral density also belongs to this universality class. 



Other arguments in favor of the unversality come from the finite volume partition 
function. Using general arguments used before in ref.u for SU(N C ), N c > 3, with 
fundamental fermions, it is possible to write down the effective partition function 
for SU(2) gauge theories with fundamental fermions and for SU(N C ) theories with 
adjoint fermions. This allows usou to derive sum rules for the inverse eigenvalues of 
the Dirac operator for an arbitrary number of flavors, and for a given value of the 
topological charge. Allresults obtained this way coincide with sum-rules derived from 
random matrix theory^. 

Questions related to the unversality of correlations in quantum spectra are dis- 
cussed in section 2. The symmetries of the QCD partition function will be analyzed 
in section 3. The main part of our program is to construct a random matrix model 
that has nothing else but these symmetries as input. This so called chiral random 
matrix model will be introduced in section 4. A summary of the analytical results will 
be presented in section 5, and numerical results for the microscopic spectral density 
in a liquid of instantons are shown in section 6. Concluding remarks are made in 
section 7. 



2. Universality of Eigenvalue Correlations in Quantum spectra 

More than three decades ago it was realized that correlations between the energy 
levels of compound nuclei show a universal behavior (see for example refHl). First, 
it was discovered that the nearest neighbor spacing distribution follows the so called 
Wigner surmise. Later it was found that not only short range correlations but also 
correlations over many level spacing do not depend on the specific dynamics of the 
system. One of the most outstanding characteristics is the stiffness or rigidity of 
the spectrum. For instance, if an energy interval contains on an average N levels, 
the variance of the distribution, £2 (AT), is not iV (as for a random sequence), but 
rather (logiV)/7r 2 . Both results can be derived from a random matrix theory with 
only the symmetries of the system as input. This provides a strong argument in 
favor of the universality of such level correlations. In particular through the work of 
Dysonta, we know that there are three distinct universality classes corresponding to 
real symmetric, complex hermitean or quaternion real hamiltonian matrices. They 
correspond to systems with time reversal invariance and integral spin or rotational 
invariance, with broken time reversal invariance, and half-integer spin systems with 
time reversal invariance and no rotational invariance. 

Much more recently, it was realized that random matrix correlations do not only 
appear in complex systems with many degrees of freedom, but also in systems with 
as few as two degrees of freedom. This was first shown by pioneering studies with 
quantum billiardsc3'EJ. Later it was found that complete chaos is not only a sufficient 
ingredient but also a necessary ingredientll. As an illustration, El we show in Fig. 1 the 
A 3 (L) statistic (an integral transform of the number variance E 2 (iV)) versus the length 
of an interval containing on an average L levels, and the nearest neighbor spacing 
spacing distribution, P(S), calculated from the (approximately) hundred lowest levels 



of the Hamiltonian 



H = i((p, - § yf + (p y - fa;) 2 ) + ^{x/arf + {y/a 2 f + g(x - yf). (9) 

For a suitable choice of the parameters a±, a 2 and the magnetic field, this Hamiltonian 
is completely chaotic for g = 0. Although in this case the time reversal symmetry is 
broken an additional symmetry allows us to choose a basis in which the Hamiltonian is 
real. As shown in Fig. 1 the spectral correlations are indeed described by the GOE. 
Only when the coupling constant is large enough to destroy all remnants of this 
symmetry, while B 7^ is kept at its original value (broken time reversal invariance), 
the spectral correlations are given by the GUE (see Fig. 2). 



Fig. 1. Numerical results for the A3 statistic (squares) and the nearest neighbor spacing distribution 
(histogram) for a nonzero magnetic field and g = 0. The full and the dashed curves show the 
corresponding results for the GUE and the GOE, respectively. 



Fig. 2. Numerical results for a nonzero magnetic field and g ^ 0. See the caption of Fig. 1 for 
further explanation. 

Two more cases of universal correlations should be mentioned. First, cross-section 
fluctuations in compound nuclei are determined by the average S matrix!. Second, 
universal conductance fluctuations: the variance of the conductivity is a given by a 
pure number times the quantum conductance 2e 2 /h. This phenomenon is closely re- 
lated to S— matrix fluctuations^, in particular to Ericson fluctuations!^ which involve 
correlations between two S matrices at different energies. 

What did we learn from the above studies? One very important implicit ingredient 
in all cases is the separation of scales: a slow varying or average scale and a rapid 



varying or microscopic scale. This allows us to express all correlations in terms of the 
slow varying scale, the average level spacing, the average S-matrix, or the average 
conductance. This would not be not possible if fluctuations on all scales were present. 
The central point is that only correlations expressed in terms of the rescaled variables 
are universal and can be described by random matrix theory. In particular, for level 
correlations it means that only correlations of the unfolded spectrum {A^} defined 
by (p is the average level density) 



are described by random matrix theory. 

Second, the microscopic correlations are completely determined by the symmetries 
of the system. In the QCD partition function, the eigenvalues show fluctuations over 
the ensemble of gauge field configurations. From what is said above we expect that 
only the microscopic correlations are universal and can be described by random matrix 
theory. In the next section we will analyze the symmetries of the QCD-Dirac operator 
and use that as an input for the random matrix model to be constructed later. 

3. Symmetries of the Dirac operator 

In this section we consider the symmetries of the QCD Dirac operator for a fixed 
external gauge field A^. The masses of the quarks are taken to be zero (chiral limit). 
In that case the Dirac operator D is chirally symmetric 



As a consequence the eigenvalues occurs in pairs ±A n , and the eigenf unctions have op- 
posite chirality. The only exception are the zero eigenvalues. For field configurations 
with topological charge v we have exactly v zero modes, all of the same chirality. 

A less obvious symmetry, which plays an important role in random matrix theory, 
is the symmetry that dictates the Wigner type of color representation of the gauge 
group. As suggested by the classical random matrix ensembles, the reality type of 
the matrix elements determines the corresponding universality class. In QCD, with 
three colors and fundamental fermions, there is no additional symmetry and the 
color representation of the gauge group is complex. Its Lagrangian is invariant under 
U(Nf) x U(Nf) (a subgroup Ua(1) is broken by the quantum fluctuations). As is well 
known this symmetry group is enlarged to U(2Nf)ac]\ for SU(2) with fundamental 
fermions, and to SU(Nf) for SU(N C ) in the adjoint representation with Nf Majorana 
fermions. This fact is closely related to the Wigner type of the representation of the 
color group. 

Indeed, for SU(2) with fundamental fermions the Dirac operator, 




(10) 



{7s,£} = 0. 



(11) 
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a 



M 2 



a 



(12) 



has an additional symmetry 



[H274r 2 color ^, J D] = 0. (13) 

Here, are the Pauli matrices acting in color space, and K is the complex conjugation 
operator. The symmetry operator satisfies the property 

(*7274T 2 col ° r ^) 2 = 1, (14) 

which allows us to choose a basis in which the Dirac operator is real. Note, however, 
that the representation of the color group is pseudo-real. Because of the symmetry 
fllcf) the symmetry group is enlarged to U(2Nf), 

For gauge theories with adjoint fermions the Dirac operator is given by 

D = i^ + if klrnllJt A™, ( 15 ) 

where the fkim are the structure constants of the gauge group. In this case the 
additional symmetry follows from the commutator 

[27274^, D] = 0. (16) 
Now, the symmetry operator satisfies 

(*7274^) 2 = -1, (17) 

which allows us to choose a basis in which the Dirac matrix can be regrouped such 
that its matrix elements are quaternion real. In this case the representation of the 
color group is real. Because adjoint fermions are formulated in terms of Majorana 
fermions, the symmetry ([IB] ) results in an SU(Nf) flavor symmetry. 

In each of the three cases discussed above the scheme of chiral symmetry breaking 
is different. For SU(N C ), N c > 3, with fundamental fermions the chiral symmetry is 
broken according to 

SU(N f ) x SU(N f ) D SU(N f ). (18) 

For SU(2) with fundamental fermions the symmetry group is enlarged to SU(2Nf) 
and spontaneously broken according to 

SU(2N f ) D Sp(2N f ). (19) 

In the case of adjoint fermions and an SU(N C ), N c > 2, color group the flavor sym- 
metry group is SU(2Nf). In this case the symmetry breaking pattern is 

SU(N f ) D 0(N f ). (20) 

The different symmetry breaking schemes are determined by the dynamics of the 
theoryB. They have to obey general constraints due to the Vafa- Witten theorem, and, 



in agreement with general QCD inequalities!!, they should give riseto pseudoscalar 
Goldstone bosons. More discussion of this point can be found in ref.H. 



4. Chiral Random Matrix Theory 

The structure of the Dirac operator in QCD is much richer than that of a Hamil- 
tonian of a completely chaotic non-relativistic quantum system. Nevertheless, it is 
possible to write down a random matrix theory that includes all symmetries discussed 
in previous section. For Nf flavors the partition function for the sector of topological 
charge v is defined by 

Z. = JVT Jd* ( $ £ ) exp(-^TYTtr, (21) 

Here, the integral is over n x m matrices T with the Haar measure. The determinant 
in this equation plays the role of the fermion determinant in QCD. We want to remark 
that the matrix 

(22) 




n ■ 



has exactly v — \n — m\ zero eigenvalues. All other eigenvalues occur in pairs ±A 
We will identify N = n + m as the volume of space time. For defmiteness we take 
m > n, and, as in QCD, we always take v<JV. The connection with the real world 
goes through the parameter E. In this model it is given by 

E = hm lim (23) 

m-fOiV-»oo TV 

According to the Banks-Casher formula@ E can be identified as the chiral condensate. 

The representation of the color group is implemented via the the matrix elements 
T and the integration measure. The number of independent variables per matrix 
element is denoted by (3 (called Dyson parameter)! 1 For 577(2) with fundamental 
fermions the color representation is pseudoreal and the Dirac operator is real. The 



integral in Eq. (|2l|) is over real matrices and we have (3 = 1. For SU(N C ) } N c > 3, 
both the color representation and the Dirac operator are complex. The matrix T in 
Eq. (|2lD is complex and (3 = 2. Finally, for SU(N C ), N c > 3, with adjoint fermions, 
the color representation is real, and the Dirac operator is quaternion real. The matrix 
elements of T are also quaternion real, and we have (3 = 4. 

An explicit realization of this partition function is for a liquid of instantons. In 
that case the Dirac operator is evaluated in the space of fermionic zero modes leading 



to the same matrix structure as in Eq. (21). However, instead of the integral over the 



*This parameter has also been introduced in the exponent, so that the average level density becomes 
independent of (3. 



matrix elements, we have an integral over the collective coordinates of the instantons. 
For a discussion and a numerical evaluation of this partition function we refer to 
setion 7. 



5. Analytical Results 

The theorical analysis of the partition function Eq. (pi]) is straightforward. The 
eigenvalue distribution is obtained by changing the integration variables according to 

T = UAV~\ (24) 

where A is a diagonal matrix with positive real matrix elements. The Jacobian of 
this transformation is0 

J(A) = U\^-^fU^ + "- 1 - (25) 

k<l k 

This Jacobian can be obtained by an explicit calculation (seeil for a discussion of the 
necessary techniques). It also follows from the fact that the Jacobian vanishes when 
two eigenvalues coincide or one eigenvalue equals zero, and that it should depend 
symmetrically on all eigenvalues. The correct powers then follow from dimensinal 
arguments. 

The integration over U and V only contributes a constant factor, so that the 
parition function is given by 



z v = c Pin jd\ x ... d \ n n \\i - ait n xt^- 1 exp (-^f ^ a 2 ,) , ^ 



where the constant Cp^ n is determined by the normalization. The joint eigenvalue 
distribution p(Ai, ■ ■ • , A n ) is just the integrand of Eq. (p6|) . We have 



Z = J dX 1 ---dX n p(X 1 ,---,X n ). (27) 

The eigenvalue distribution follows by integrating over all eigenvalues except one 

p(Ai) = J dX 2 ---dX n p(X l ,---,X n ). (28) 

Integrals of this type have been widely studied in the context of random ma- 
trix theory. One very powerful method is the orthogonal polynomial method Dyson, 
Wigner and Mehta (see the book of MehtaEa for references). This method is partic- 
ularly suited for (3 = 2, when the integrals follow immediatly from the orthogonality 
relations of the orthogonal polynomials. For (3 = 1 and (3 = 4 it is still possible to 
do the integrals, but the resulting expressions are much more complicated. Instead 
of classical orthogonal polynomials, one has to invoke skew orthogonal polynomials. 
The general method was developed by Dysonll and by Mahoux and MehtaH and 



applied to many different ensembles by Nagao, Slevin and WadatS Unfortunately, 
the generalized Laguerre ensemble for (3 = 1 was not analyzed before, so we perfomed 
the analysis in ref.e3. The symplectic generalized Laguerre ensemble has not yet been 
studied, but using similar methods it is possible to obtain an explicit expression for 
the level density. 

First we discuss the simplest case, (3 = 2. This case is also known as the unitary 
generalized Laguerre ensemble and was studied inBac3c3E!l. The latter authors also 
evaluated the microscopic limit. This ensemble was put in the context of QCD in 
refHI, where the microscopic spectral density for v = was derived. The microscopic 
spectral density for arbitrary v was first given in refc3. 

The result for the microscopic spectral density @ for the chGUE {(3 = 2) 



Ps( z ) — —^-(Jnj+uC^z) - JNf+v+lQ^z) JNf+v-li^z)). 



(29) 



It depends only on the combination Nf + v which agrees with the sum rules derived 
by Leutwyler and Smilgap. It is also possible to obtain explicit expressions for the 
two level correlation functions. For results we refer to the literature^. 

The expressions for SU(2) with fundamental fermions (J3 = _L) are much more 
complicated. In this case we find the microscopic spectral density^ 



Ps{z) 



y y r°° 

-J 2a+ i(zy) + -/ dw(zw) 2a+1 e(z 
4 2 Jo 

wJ 2a (zT,)J 2a _ 1 
x — - 



(1 d 
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WS) - zJ 2 a- 1 (zy)J2g 

zw) 2a (z 2 
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where a is the combination 



N f 



v 
2' 



(30) 



(31) 



Eqs. ( |29| ) and (|30| ) can be used to derive Leutwyler-Smilga sum rules. However, 

type of integrals 
for the quantities 



there is a more direct way to 
follow from the Selberg formula 



fin these results. Exactly, thes 
. Sum rules have been obtained^ 



\n\i=n%i^--i=-n v ni 



XI--- X 2 



(32) 



and 



S, 



pq \ ^ A 4 • • • A 4 A 2 • • • A 2 

\ all Aj different mi "my'^i y n q 



(33) 



The integrals can be calculated by deriving a recursion relation. It is also possible 
to derive recursion relations for more complicated sum rules. However, then the 
recursion branches off and cannot be solved analytically any longer. 



The results for the simplest family of sum rules are given by 



Sp = > (34) 

which, in the limit N — > oo, simplifies to 

(jv 2 s^ iy+ 1 + 2(^-1)/^) 

p 2>! T(v + p+l + 2(N f -l)/p)' K ' 

This result for (3 = 2, and for (3 = 4, p = 1 and Afy = 1 or iVy = 2. were derived 
previously with the help of the finite volume static partition functional. Results for 
all other case, in particular for SU(2) with fundamental fermions, were first derived 
with the help of random matrix theory^. Presently, we have shown that the p == 1 
sum rules can also be obtained starting from the finite volume partition function^. 

Note, that for Nf = 1 the parameter (3 drops from the equation. Indeed, in this 
case the finite volume partition function is solely determined by the Ua(1) anomaly 
and is not dependent on the Wigner type of the color representation. There are no 
additional degrees of freedom that, for more than one flavor, contain information 
on the color representation. This implies that the microscopic limit of the spectral 
density is not determined uniquely by the finite volume partition function. 

The results for S pq can be derived in a similar fashion. We only quote the resu 

r(a + |)r(a + p)T(a + | +p + n) , 2 2 y+2 P 
Spq = Npq T(a + l + q + 2p)T{a)T{a + l + n) \ n S ) ' (36) 

where the factor N pq is defined by 



N pq = T7? : rr (37) 

p\q\[n — p — q)\ 

and 

Q= , +1+ 2 w- 2) - m 

P 

These sum rules depend on (3. Indeed, the effective low energy field theories are 
different for different values of (3. Even for Nf = 1 they depend on (3 (the integrals 
should be convergent). Since such sum rules cannot be derived from the finite volume 
partition function these results are completely consistent. 

The result (^) was derived from the finite volume partition functional only for 
(3 = 2, p = 1 and q = 0. 



6. Numerical Results 



It would be very interesting to check the results of section 5 by lattice QCD 
simulations. However, in practise, it is not possible to obtain a sufficient number of 



eigenvalues of the Dirac operator close to the chiral limit (see however ref.0, where the 
spectrum near zero virtuality has been obtained for quark masses that are much larger 
than the smallest eigenvalues, and in the present context, corresponds to Nf = 0). 

Now, remember that the correlations discussed above are based a random Dirac 
operator. When we consider the fluctuations of the eigenvalues of the Dirac operator 
over the ensemble of the gauge field configurations, it seems clear that if the ran- 
dom matrix results are reproduced by cooled configurations, they are certainly found 
for the original uncooled configurations. The sum over all cooled configurations is 
well known: it is a dilute gas or liquid of instantons. Although, it does not have 
confinement, the bulk of hadronic physics can be reproduced bv this modefi 

The partition function for a liquid of instantons is given bycj'19 

z = / n^n det ( TtT + m /) ex p(- 5 giue), (39) 
i=i f 

where the integral is over the collective coordinates of each instanton, and the Dirac 
opertor is calculated in the space of fermionic zero modes with matrix elements de- 
noted by Tij. The gluonic action, S s i ne , contains the interaction between the instan- 
tons approximated by a sum over all pairs. The total number of instantons in N. 
In general, the number of instantons is not equal to the number of anti-instantons. 
However, in the numerical calculations to be presented below, this will be the case. 
Further details of this model and the calculations presented in this section can be 
found in ref.il. 



Fig. 3. A histogram of the total spectral density of the Dirac operator in a liquid of instantons 

In the partition function fl39|) it is straightforward to change the number of colors 
from N c = 2 to N c = 3. In fact, we use the same computer program that covers both 
cases. At present it is not yet clear how to write down an instanton partition function 
for adjoint fermions. For a discussion of the difficulties related to this issue we refer 
to refH 

The numerical results have been obtained for 32 instantons and 32 anti-instantons 
confined to a box of 2.37 3 x 4.74A^ D . For N f = 0, the average was over 100,000 
configurations. For all other flavors the average was over only 20,000 (correlated) 
configurations. A histogram of the results for the total spectral density is shown in 
Fig. 3. The results are for two and three colors and one, two and three flavors (see 
label of the figure). In Fig. 4 the region of the spectrum near zero virtuality is put 
under the microscope. We show a histogram of the unfolded spectrum in terms of the 
microscopic variable z, so that the average spacing between the eigenvalues is one. 



Fig. 4. The microscopic spectral density. The full line reprsents the numerical results. Random ma- 
trix result of the chGOE and the chGUE is shown by the dotted and the dashed curves, respectively. 

In the same figures we also show the random matrix results for the chGOE (dotted 
curve) and the chGUE (dashed curve). The only deviation is that the oscillations in 
the chGUE random matrix results seem not to be completely reproduced. At the 
moment it is not clear whether that this is a real effect, of that it is a consequence of 
the finite size of the system. 



7. Conclusions 



We have argued is that the average miscroscopic spectral density of the QCD 
Dirac operator follows from its symmetries and a maximum entropy principle. The 
latter principle uniquely determines the random matrix theory obeying all constraints 
imposed by the symmetries. At present, it is a conjecture that the the fluctuations 
of the eigenvalues near zero virtuality over the ensemble of gauge field configurations 
as weighted by the QCD partition function belong the the same universality class as 
the corresponding chiral random matrix theory. We have presented several strong 
arguments in favor of this conjecture. First, results from the Hofstadter model show 
that the random matrix model has a wide basin of attraction. Second, the triality 
of the classical random matrix ensemble also occurs in the chiral random matrix en- 
sembles and coincides with the different schemes of chiral symmetry breaking. Third, 
numerical calculations for instantons agree with the random matrix result. Fourth, 
all spectral sum rules that can be obtained from the finite volume static partition 
function, which coincides with the full QCD partition function in its region of appli- 
cability, are also given by the chiral random matrix theory. 

The finite volume static partition function alone does not determine the micro- 
scopic spectral density uniquely. The most striking example is that of one flavor. In 
that case the effective theory does not depend on the color representation, whereas 
the microscopic spectral density is different in each of the three universality classes. 

To address the question which effective field theory contains the information of 
the full microscopic spectral density, remember that the resolvent is defined by 

G M?^>- (4o) 

It can be expressed as follows in terms of the partition function 

/ 1UA, ■■■-) A 

dz 



W--) = 4z ( %»l*l m% + ™f))- (41) 



\lKA7, • / 

In ( flop and (|41|) the sum and product is over the positive eigenvalues only. Therefore, 
the spectral density for Nf flavors can be obtained from a partition function of Nf + 1 
fermionic flavors and one complex scalar flavor. 
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